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Abstract: Lop Sea, located at the east end of the Tarim Basin, Northwest China, dried up permanently, 
which is the terminal lake of the Tarim River. Lop Sea was considered as the lake basin of Lop Nor since 
Quaternary. However, the possibility that Lop Nor was away from the Lop Sea in historical time is crucial 
to be discussed to interpret the proxy records in sediment profiles. To obtain a general view of the evolution 
of Lop Nor and Lop Sea in a historical period, several approaches were adopted in this paper. First, the 
Qianlong Thirteen-Row Atlas, an ancient imperial atlas of the Qing Dynasty, which was completed around 
1760, indicated that the Tarim River formed a relatively large lake at its modern upstream region. Second, a 
Digital Elevation Model (DEM) with a 10-m spatial resolution and a relative precision of 0.42 m was derived 
from TanDEM-X/TerraSAR-X satellite image pairs using the interferometry method, which was verified 
using ICESat-GLAS laser footprints and a local DEM acquired by a drone. Finally, based on the spatial 
analysis of historical documents, expedition reports, sediment profiles and archaeological evidence, it can 
be deduced that the lacustrine deposition was discontinued in the Lop Sea. Six episodes in the evolutionary 
history of the drainage system in eastern Tarim Basin were summarized. The proved depositional condition 
variations could be used for future interpretation of proxy records in sediment. The high-accurate DEM 
provided a reference for the location of further fieldwork in the Lop Sea. The method proposed in this 
paper may be efficient for the research of inland lakes or rivers in global arid regions. 
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1 Introduction 


Lop Nor, located at the east end of Tarim Basin is dried up permanently because of climate 
change. At the juncture of the westerly and the monsoon (Chen et al., 2008), Lop Nor is 
important for studying the impacts of regional paleoclimate changes on human activities. The 
location of the terminal lake of the Tarim River, assigned as Lop Nor, had been controversial 
over the last century. The name of Lop Nor initially arose on Qianlong Thirteen-Row Atlas of the 
Qing Dynasty (Wang and Liu, 2007), which was designated as Qing Lop Nor (Fig. 1) in this 
paper. From 1876 to 1916, instead of the location illustrated on the map, the terminal lake of 
the Tarim River was found at Kara-koshun (Fig. 1) by various famous expeditions led by 
Hedin (1905), Huntington (1907a), 
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Przhevalsky (1999), and Stein (1921). The debate about the location was mainly between 
Przhevalsky (1999) and Hedin (1905). Przhevalsky (1999) thought Kara-koshun was the Qing Lop 
Nor due to immature cartography. Hedin (1905) insisted that Qing Lop Nor should be at the eastern 
end of a wadi Kuruk-daria (Fig. 1). Around 1921, instead of going southward, the Tarim River 
flowed to the Kuruk-daria at Temenpu (Fig. 2) and reached the northwest part of the Lop Sea (Lop 
Nor 1931 in Fig. 1) (Hedin, 1905; Hörner and Chen, 1935; Chen, 1936;) supporting Hedin's 
judgement. 
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Fig. 1 Sketch map of Lop Region. Lakes and courses according to Hedin (1905), Hörner and Chen (1935), and 
Qianlong Thirteen-Row Atlas (Wang and Liu, 2007). The names of the ancient city ruins follow Stein (1921). The 
term 'Lop Sea' follows Hörner and Chen (1935) to distinguish Lop Nor in other periods. Contours are derived from 
SRTM (Shuttle Radar Topography Mission, http://www2.jpl.nasa.gov/srtm/). 


Lake Lop Nor might be completely dried up after 1973, based on the interpretation with decades 
of remote sensing data (Xia et al., 2007). With more geo-informational data and field materials, the 
Lop Sea was thought to have been the permanent terminal lake of the Tarim River (Zhou, 1978; 
Yuan and Yuan, 1998; Xia et al., 2007; Dong et al., 2012). Under the continuous lacustrine deposit 
assumption, sediments in the Lop Sea were sampled and interpreted (Ma et al., 2008; Jia et al., 
2011; Zhang et al., 2012; Liu et al., 2016) to infer the regional paleoenvironment. However, no 
convincing evidences were found to validate this assumption. On the contrary, Chinese literatures 
implied that the courses and the terminal lake of the Tarim River changed several times during the 
historical period (Xi, 1992). Lacustrine deposition discontinuity is often discovered and defined by 
dating materials at high sample resolution (Coianiz et al., 2015; Ivory and Russell, 2018) or 
combined with sedimentary facies analysis (Jiang et al., 2015). The available dating materials are 
insufficient at a suitable sample resolution to identify short lacustrine discontinuity. With the 
development of remote sensing, the formation of the elevation of Lop Nor has been demonstrated 
to be tied to the composition variation of the salt crusts (Li et al., 2008). Digital Elevation Model 
(DEM) can provide geomorphological evidence for the reconstruction of lake evolution (Reinhardt 
et al., 2008). 

In this paper, a spatial analysis of the published researches was carried out to deduce the 
evolution of the Lop Sea during the historical time. Qianlong Thirteen-Row Atlas (Wang and Liu, 
2007) and other historical documents was used to verify the location of Qing Lop Nor. Besides, a 
highly precise DEM was produced to describe the topography of the Lop Sea and surrounding 
regions. Based on the interpretation of DEM combined with expedition reports, sediment profiles, 
and involved radiocarbon dating ages, possible lacustrine deposition discontinuity of Lop Sea are 
discussed, and six typical episodes of the drainage system evolution process are summarized. 
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2 Study area 


The Lop region studied in this paper (39°-41°N, 89°—92°E) is the lowest section of the Tarim Basin 
with an average elevation of 790 m. North of the Lop region is Kuruktag Mountains, south is 
Astintag mountains, and west is Taklimakan Desert (Fig. 2). The western desert of the Lop region 
is called the Lop Desert, while the eastern salt crust is the Lop Sea. Based on the current 
hydrological settings, after being fed by rivers from Khotan, Yarkand, Shache, Kashgar, Aksu, and 
Kuche, the Tarim River flows south from Arghan and runs into Kara-buran (Fig. 2). The runoff 
volume of the Tarim River varies with seasons (Zhang et al., 2008) and years (Xu et al., 2005). As 
a result, the water coverage area of the lower reaches of Tarim River varies greatly (Abdimijit et 
al., 2016). Another northern tributary, the Kongque-daria River, terminates at Ak-supu (Fig. 2). The 
Cherchen-daria River, which is another supplement of Kara-buran and Kara-koshun, did not feed 
the terminal after the 1980s due to increasing demand for irrigation. The region is uninhabited now, 
except for a potash fertilizer mine, lying in the center of the Lop Sea. 

The annual average temperature is 11.6°C, with the maximum summer temperature higher than 
40.0°C, and the minimum winter temperature below —20.0°C. The average annual precipitation is 
less than 30 mm, while the potential annual evaporation is 3500 mm (Mischeke et al., 2020). The 
prevailing northeastern wind has shaped the Yardangs landform around the Lop Sea. 
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Fig. 2 Drainage system of the Tarim Basin at present. Tinted SRTM DEM (Digital Elevation Model) map 
represents the topography of the Tarim Basin, see the color bar for detail. 


3 Materials and methods 


The famous Wu-Chang-Fu Map (Han and Lv, 2006), which guided the explorers a century ago, was 
adapted from a private imperial map named the Qianlong Imperial Atlas, or the Qianlong Thirteen- 
row Atlas (Lu, 2008). The emperor of Qing, Qianlong, sent his subjects led by HE Guozong and 
Minggatu, along with two Portuguese Jesuit missionaries, José de Espinha and Felix da Rocha, to 
survey the astronomical coordinates of the main settlements of Xinjiang in 1756 and 1760 (Guo 
and Li, 2003, 2005). Based on the survey results, the cartography work was completed by a French 
missionary, Michel Benoist (Wang and Liu, 2007). The prime meridian was the one going through 
the Forbidden City, Beijing. Just as its name implies, the sheet line system with lines of integral 
latitude and longitude divided the atlas into thirteen rows according to latitude. The eighth row, 
western #3 of the Qianlong Thirteen-row Atlas contains the geographical elements in this paper. 
The mother sets of its copper plate, which were engraved by 1775, were found at the Palace 
Museum in 1925 and were republished in 1931 (Wang and Liu, 2007). The republished map was 
scanned by a Sony a7RM2 camera (Sony Corporation, Japan) with a Sigma 85 mm F1.4 lens 
(Sigma Globle Vision, Japan), and geocoded by the ArcGIS tool Georeferencing (Environmental 
Systems Research Institute, Unite States) using the ‘project transform’ method, taking the 
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intersections of latitude and longitude as control points. 

Then, the maps were validated with historical literature. The latitudes and longitudes of 46 main 
towns of Xinjiang and central Asia were recorded as the control points for cartography (Guo and 
Li, 2003, 2005). Here, six of the towns near Lop Nor, Luntai, Kuche, Luke-qin, Turpan, Hara-shar, 
and Korla (Table1; Fig. 2), which were in the interesting subdivision map, were selected to evaluate 
errors. The current geographical coordinates of these towns were extracted from Google Earth. 

The public global DEM data sets, such as SRTM and ASTER (Advanced Spaceborne Thermal 
Emission and Reflection Radiometer) were not good enough to interpret the drainage system. Thus, 
the Synthetic Aperture Radar (SAR) data and TanDEM-X/TerraSAR-X Coregistered Single look 
Slant range Complex (CoSSC) were used in this research, which was supported by the German 
Aerospace CenterTanDEM-X Science Coordination, with the project names of OTHER6906 and 
ATL HYDR7333. To cover the region of interest, 28 pairs of images acquired in 2015 were selected 
(Fig. 3). The synthetic aperture radar interferometry (InSAR) DEM Workflow module of 
SARscape 5.2 (SARmap company, Switzerland) was adopted for data processing, with SRTM V4 
as an initial reference DEM (Sun et al., 2016). The output spatial resolution was set to 10 m. 
Additional ground control points (GCPs) were extracted in these completed DEMs for other image 
pairs. The final product was mosaicked by the ENVI (Environment for Visualizing Images) 
softwareseamless mosaic module without feathering or averaging. 

To control the global precision, screened ICESat/GLA14 data from 2003-2009 (Li et al., 2017) 
is chosen as the GCPs because its height accuracy is usually better than 1 m (Gonzalez et al., 2010). 
The coordinates were transformed from Topex/Poseidon to WGS84 (World Geodetic System 1984; 
Li et al., 2017). The gravity anomaly was amended by the Earth Gravitational Model EGM2008 
(Pavlis et al., 2008, 2012). Due to the valid footprints of GLAS (Geoscience Laser Altimeter 
System) spread sparsely in the study region, 8 pairs of CoSSC, where footprints were in enough 
quantity and extensive distribution, were selected as the starting control images (Fig. 3). 
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Fig.3 Distribution of Co-registered Single look Slant range Complex (CoSSC) image pairs and Geoscience Laser 
Altimeter System (GLAS) footprints, with SRTM as the background elevation reference 


The DEM accuracy was assessed from two points of view, including the absolute root mean 
square error (RMSE) and relative local error. To estimate the global RMSE, all the valid GLAS 
footprints were referred as the truth elevation value of the ground. To evaluate the local 
performance, the ruined fort LK (89°40'4.75"E, 40°05'39.40"N; Fig. 1) is selected as a validation 
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site. In 2016, aerial photogrammetry Structure from Motion (SfM) (Mancini et al., 2013; Lucieer 
et al., 2014a, b; Wallace et al., 2016) fieldwork was conducted to obtain a 3D model of the ancient 
fort. Photos were taken by DJI Phantom 3 Pro (Shenzhen DJI Sciences and Technologies Ltd., 
China) 30 meters above land, with more than 80% overlap of an approximately 200 mx200 m area 
covering the ruin. 


4 Results 


4.1 Validation of the digitalized Qianlong Thirteen-Row Atlas 


In 260 years ago, the vertical accuracy of documented astronomical coordinates (Table 1) was quite 
high. Errors in longitude are greater than those in latitude, which is a normal regulation for 
astronomical surveying. The worst bias in longitude occurred in Hara-shar. 

The error might also come from cartography and digitizing. The forward RMSE of geocoding is 
0.01556°, which is equal to 0.93’, with 9 GCPs spanning the map. Comparing coordinates extracted 
from the geocoded map and the current map (Table 1), it is found that there are significant errors at 
Luntai and Kuche, which is impossible in digitizing, because the longitude bias is close to one degree. 
The only explanation is that Michel Benoist, the author of the map, made some adjustments to 
integrate the surveying results, but we don't know how these errors occurred. Fortunately, the 
coordinates of Hara-shar and Korla are consistent with their current geolocations with an acceptable 
bias. 


Table 1 Coordinates of main towns near Lop Nor 
Town Latl Lonl Lat2 Lon2 Lat0 Lon0 Latl-Lat0 Latl-Lat0 Lat2-Lat0 Lat2-Lat0 


Luntai 41°44'00" 84°16'27.00"41°32'31.20"83°27'22.80"41°46'30.00"84°15'15.00" —2'31.80" 1112.00" -14'0.60" —47'52.20" 
Kuche 41°37'00" 82°54'27.00" 41°24'9.00" 82°00'37.80"41°42'57.60"82°57'1.20" —5'57.60"  —5'56.40" —18'48.60" -56'23.40" 
Luke-qin 42°48'00" 90°12'27.00"42°46'15.60'90°04'13.80"42°44'41.4089°45'34.80"  3'18.60"” 26'52.80" 134.20" —18’39.60" 
Turpan 43°04'00” 89°38'27.00"42°56'34.20"89°37'32.40"42°56'59.40"89°11'19.20"  7'0.60" 278.40" = -0'25.20" = 26'13.20" 
Hara-shar 42°07'00" 87°06'27.00'42°11'39.00'86°31'43.20"42°03'46.8086°34'10.80" 313.20" 3216.20" 7'52.80" —2'27.00" 
Korla 41°46'00" 86°27'27.00'41°30'28.8085°51'45.00” 41°44'3.60” 86°10'21.60"  1'56.40” 17'5.40" —13'34.80" —18'36.60" 


Note: Lat, latitude; Lon, longitude. Lat! and Lon! are values recorded in documents, where longitude, subtracted by the longitude of the 
Forbidden City 116.390888°E, has been transformed as the Greenwich Prime Meridian. Lat2 and Lon2 are extracted from geocoded 
geotiff file of Qianlong Thirteen-Row Atlas. Lat0 and Lon0 are the current geo-coordinates extracted from Google Earth. 


4.2 DEM and accuracy evaluation 


The distribution of error was projected on the tinting DEM to show the spatial features (Fig. 4). 
The histogram (Fig. 4) shows that the errors nearly yield a normal distribution. The global RMSE 
of the DEM is 0.42 m. Errors in the eastern playa are minute, usually less than 1 m, while the errors 
in the desert and Yardangs, are sometimes substantial. The possible reason is that the geolocation 
of data is not rigorously aligned, so the elevation difference will be dramatic if the local landscape 
varies sharply, such as it does in the Yardangs region. Overall, the DEM is accurate enough to tint 
with a 1-m hypsometric interval. 

The InSAR DEM (Fig. 5a) was validate with the DEM (Fig. 5b) in a 2-cm spatial cell with a 
relative vertical accuracy of one decimeter derived from the SfM. Five spots in LK (Fig. 5) were 
selected for quantitative comparisons, and the validation results are summarized in Table 2. Since 
they are in different elevation systems, it makes sense to compare the relative height rather than the 
absolute height. The flat ground to the east of the ruined fort (Fig. 5) is taken as the local base 
height. 

Although the InSAR DEM roughly presented the topography features of the ruined ramparts, the 
relative height of three points on the tops of the ramparts, which are higher than the assigned base 
height, are all underestimated, up to 2.8 m at the southwest wall of the ruin LK. There is a horizontal 
deviation of 20 m between two images at the southeast wall of the ruin LK, resulting from the 
resolution differences. The elevation of the two depressions inside the ruin is also underrated. The 
relative vertical error in such complicated terrain could be evaluated as 3 m, as shown in Table 2. 
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It is not difficult to understand these differences. On the one hand, the 10-m spatial resolution limits 
the validity. On the other hand, both walls and depressions, are isolated islets with a sharp height 
variation on their margin, where the interferometric fringes would discontinue easily on an 
interferogram (Fritz et al., 2012). The errors are due to the limitation of the radar system. 
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Fig.5 Comparison of DEM for the ruin LK from different resources. (a) nSAR DEM; (b) aerial photogrammetry 
SfM DEM. The color ramp is the same as in Figure 4. The both are projected on WGS84 UTM Zone 45N. NE, 
northeast; SW, southwest; SE, southeast. 


5 Discussion 


5.1 Interpretation of nSAR DEM 


The figure about the area and boundaries of Lop Sea has not been ever provided. The spatial 
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Table 2 Relative elevation differences of the two DEMs 


Elevation (m) 


Site Error 
InSAR SfM InSAR relative SfM relative 

Reference zero height 793.600 791.739 0.000 0.000 - 
Northeast wall 800.000 799.177 6.400 7.438 —1.038 
Southwest wall 799.800 800.775 6.200 9.036 —2.836 
Southeast wall 797.400 795.837 3.800 4.098 —0.298 
Northern depression 793.600 790.307 0.000 —1.432 1.432 
Southern depression 795.800 791.679 2.200 —0.060 2.260 


Note: InSAR, synthetic aperture radar interferometry; SfM, Structure from Motion. 


difference between the west shore of Lop Nor 1931 (Hörner and Chen, 1935) and the 790.0-m 
contour in InSAR DEM is approximately 5’ in longitude, and negligible in latitude. It is reasonable 
for astronomical surveying results to contain a larger error in longitude than in latitude, although 
they can satisfy the precision requirement. Thus, the boundaries of Lop Sea with 790.0 m could be 
confirmed. To validate the remaining bank at the same elevation of approximately 790.0 m, three 
field sites were checked to the north, east, and south (Shao et al., 2011). The bank is confirmed on 
sites, where the inward region of the salt crust is mixed with fine sand, the outer region is clay bluff 
to the east, and the beach covered with gravel to the north. Meanwhile, the south site is covered by 
an alluvial fan, suggesting that the ancient southern shoreline may have been further south than 
nowadays. 

On the one hand, the water level could not be higher than the elevation of the ruins, LE (790.4 
m), LA (792.8 m), and LK (793.8 m) in the Loulan period (Fig. 6) (approximately 25-400 AD (Lv 
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Fig. 6 Interpretation of InSAR DEM. It is projected on WGS84 UTM Zone 46N. The color ramp is the same as 
in Figure 4. Pit I is the location of Liu et al. (2016). Pit II is the location of Ma et al. (2008). The red dotted line 
square demonstrates the spatial range of Loulan. 
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et al., 2010; Xu et al., 2017)). Thus, the maximum water level of the Lop Nor in the Loulan period 
was about 790.0 m. The area of the Lop Nor in the Loulan period could be extracted from the DEM, 
which is 8156.0 km? with a water volume of 16.50 km?. This may be the region of the Puchang Sea 
near Loulan city recorded in the Book of Han, which comprises an area of 126 km by 126 km in 
longitudinal and latitudinal directions; however, this is an estimated data and measurement method 
is unknown. On the other hand, the salt crust beyond the west strand (Fig. 6) (Shao et al., 2012) 
with an average elevation 791.4 m implies that the Lop Sea at a previous time might have been 
larger than the one enclosed by the 790.0-m contour. 

A distinct broad channel 2 (CN2) connects the Lop Sea and Kara-koshun (Fig. 6). Based on the 
result of InSAR DEM, the width and depth of CN1 (the wide southeastern channe 1) and CN2 can 
be obtained. The widest part of CN2 is approximately 2.9 km at its beginning in the northeastern 
Kara-koshun, while the narrowest part is approximately 0.4 km in the middle. The average depth 
of CN2 is over 1 m. In the east of Kara-koshun, beside CN2, there are some fork-shaped trenches 
stabbing toward Lop Sea. The northern strand of Kara-koshun is clear and continuous. The 
elevation in the west and south varies smoothly, where no distinct channel is identified. 

North of the Kara-koshun depression are the other two depressions, the Loulan depression and 
the central depression, which were first discovered and named by Hedin in his leveling survey from 
LA to Kara-koshun in 1901 (Hedin, 1940). They are now clearly identified in a topographic map 
for the first time. 

Two broad and deep-undercutting channels CN1 and CN2 to the south of LA, SR1 (Southern 
River) and SR2 (Southern River), extend into the Loulan depression (Fig. 6), suggesting that this 
depression accepted intermittent flooding after its desiccation. The continuous eastern levee of the 
Loulan depression and CN1 connecting the Lop Sea (Fig. 6) could be used as evidence for this 
inference. The average elevation of the levee is 790.6 m, while the lowest elevation of CN1 is 788.6 
m, which indicates that the depression was discharged into the Lop Sea through CN1 when the 
water level reached 789.0 m, and the maximum water level in Loulan depression was approximately 
790.0 m. In terms of the water level at 790.0 m, the possible lake area of the Loulan depression 
would be 354.3 km’ with a volume of 0.41 km’. 

The lakes in arid areas are sensitive to climate change and human activities, and the changes in 
lake historical water level affect human activities. Under the global warming environment, lakes in 
arid areas began to shrink, exposing the original lake basin. The topographic survey method used 
in this paper can be adopted to study ancient lakes in arid areas. 


5.2 Location of Qing Lop Nor 


The historical documents showed that no astronomical surveying was conducted in the south or 
east of Lop Nor, which means that the measured lakes are the consequence of traverse surveying 
from Korla to the northwest. Due to an unknown surveying method for the lakes, a reasonable 
assumption should be made to estimate the bias. Regardless of the adopted methods, the location 
must be derived from measured angles and distances. Considering the surveying task of Espinha in 
1759, including the longitude from 86°E (Korla, China) to 69°E (Tashkent, Uzbekistan), in just one 
year, it was difficult to guarantee the precision of the details. The errors of angles and distances are 
supposed to be 10° and 10% of the total distance from the origin, respectively. The error ellipses of 
the river mouth (Fig. 7) were drawn to represent the influence of these errors on a survey. 

On the Qing map, the Tarim River converged with Kongque-daria at 90 km south from Korla 
and flowed 110 km to feed Lop Nor (Fig. 7). It was reported that Lop Nur change into a group of 
small lakes during the Qing Dynasty (Yang et al., 2006). The error should not be too large, just as 
the demonstrated error ellipse. This location is in well agreement with the records by a Chinese 
geographer, XU Song, in his book Hydrology of Western China, in which Lop Nor was reported 
287 km away from Hara-shar. If this distance was counted as the route along Kongque-daria from 
Hara-shar, the water coverage area of Lop Nor in the Qing map should be about 2536.0 km?. 

There was no river to connect the surrounding small lakes and Lop Nor, perhaps they were fed 
by seasonal floods. There was a lake, Ghas Nor, to the south of Lop Nor. This lake might be Kara- 
buran, according to the error ellipse, which agreed with Huntington (1907b). Between Ghas Nor 
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and Lop Nor, no channel directly connected these two lakes, except for two courses named Kubu 
and Nu-qi-tu-da-ba-da-han (Fig. 7). The latter might be the rudiment of the group of lakes, Avullu- 
kol, Kara-kol, Tayek-kol, and Arka-kol, which were thought to be the remnant lakes of Lop Nor in 
the Qing map by Hedin. (Figs. 1 and 7). 

Hedin (1905) thought that Qing Lop Nor was located in the Loulan depression, but with the 
InSAR DEM, it can be seen clearly that the shape of the Loulan depression was different from that 
of Qing Lop Nor (Fig. 7). In addition, its area (354.3 km?) was not large enough to hold such a 
huge lake (about 2536.0 km”). No suitable depression agreed with the shape of Qing Lop Nor along 
Kuruk-daria (Fig. 7), which was against the conjecture of Huntington (1907b). Their views about 
the location of Qing Lop Nor could not be supported by topographical evidence. 

Qing Lop Nor should be on the place illustrated in the ancient map. Besides, the water coverage 
area of Qing Lop Nor (about 2536.0 km?) is close to that of Lop Nor 1931 (about 2445.0 km’), it 
is highly possible that Qing Lop Nor was the terminal lake of the Tarim River in the Qing Dynasty, 
indicating that Lop Sea might get smaller or ephemeral water supply during this period (Mischke 
et al., 2020). 
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Fig. 7 Error analysis of Qing Lop Nor projected on WGS84 UTM Zone 45N. Error ellipses of river mouth are 
drawn by red lines. Map of Hedin 1896 is compared, with SRTM as an elevation reference background. 


5.3 Spatial evidence for the CN2 formation dating 


Another debatable issue in Lop Nor evolution research is the active duration of CN2 (Fig. 6). It 
was thought to be formed around 1500 a ago (Yuan and Yuan, 1998), while this conclusion is 
challenged by the current topographical evidence. Kara-koshun was reported to overflow 
northwards to form a new marsh Yangi-kol (Fig. 1) in 1901 (Hedin, 1905). If CN2 with an elevation 
of 788.5 m (Fig. 6) was formed before 1901, Kara-koshun could not maintain the water level at 
790.5 m to overflow in dry winter. In addition, the leader of Abdal (Fig. 1), declared in 1890 that 
no channel could be found in the northern Kara-koshun after he travelled around the marsh seeking 
better homesteads for his people (Pjevtsoff et al., 2013). However, Kara-koshun discharge to the 
Lop Sea by this channel, which can be confirmed by Huntington's travel notes when he traversed 
the playa from the south to north at 90°48’E (the route of Huntington in January 1906; Figs. 1 and 
6). The observed ‘whiter, fresher and slightly damp’ salt crust (Huntington, 1907b) demonstrated 
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that this region should have been inundated not long before. Later the channel was reported deserted 
in 1931 (Chen, 1936; Hörner and Chen, 1935). Thus CN2, which was the only drainage channel 
from Kara-koshun to the Lop Sea, might be formed between 1901 and 1906. The corresponding 
terminal lake area in 1906 (inferred water area in 1906; Fig. 6) could be inferred according to the 
expedition route, covering about 1580.0 km’. The evidences also indicated that the Lop Sea might 
be subaerial before 1901. 


5.4 Sediment evidence for the discontinuous lacustrine deposition 


Among the published sediment researches, Pit I (Zhang et al., 2012; Liu et al., 2016) and Pit II (Ma 
et al., 2008) represented good grain size sequences (Fig. 8) at their locations (Fig. 6). In Pit I, the 
OSL (optically stimulated luminescence) ages at depths 37 and 75 cm were 500 (+100) and 1800 
(+200) a BP, respectively (Zhang et al., 2012). A dramatic grain size peak appears at a depth of 
around 46-54 cm. While in Pit II, the sediments at a depth of 62—80 cm are dark gray silty clay 
with plant oddments, whose radiocarbon ages were 871 (+45) (at 70 cm) and 800 (+40) a BP (at 75 
cm). The grain size increases rapidly in the overlying above 28 cm. Neither the depth-age relations, 
nor the grain size patterns, are consistent in these two profiles, which is hard to explain under the 
stable lacustrine deposit assumption. 
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Fig. 8 (a) Grain size at Pit I (Liu et al., 2016; Zhang et al., 2012) and (b) Pit II (Ma et al., 2008), re-interpreted in 
this paper. The dating of Pit I was OSL from Zhang et al. (2012). Dating of Pit II was radiocarbon age of reed root 
in sediment. The blue rectangle represents the period with water supply. 


Grain size can reflect the change of river mouth in Lop Sea. In 1896 and 1900, Kuruk-daria (Fig. 
1) was reported as a dry channel without any water supply (Hedin, 1905, 1940). However, the river 
mouth was discovered at the end of Kuruk-daria, near Pit II in 1931 (Hörner and Chen, 1935) and 
1934 (Chen, 1936; Hedin, 1940). Besides, near the river mouth, a layer of soil sediment at the 
bottom of the lake, whose underlying was the solid salt crust, was measured about 25 cm (Chen, 
1936), which coincides with the situation in Pit II, where the coarsest grain occurs at 0-28 cm depth. 
It was believed that this layer of soil was the result of the latest water supply after 1921 (Chen, 
1936). This modern case presented an example for the increasing grain size with the changing of 
the river mouth and the discontinuous lacustrine deposit in Pit II. 

An identical event might be recorded at a depth of 46-54 cm in Pit I, where the grain size 
increases rapidly. Pit I is located in the alluvial fan of CN1. It is consistent with the wet season 
indicated by the dating results of the samples from the Loulan area (Li et al., 2018). Li et al. (2018) 
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collected natural plant relicts and samples from active traces of human in the surrounding area of 
Loulan ancient city for AMS (accelerator mass spectrometry) 'C chronological analysis. During 
1260-1450, the Tarim River continued to drain into the Lop Sea through the Loulan River. 

According to the spatial reference supported by the reliable DEM in this study, if the Tarim River 
flowed through SR1 and SR2 into the Loulan depression, it would be highly possible that it was 
drained into Lop Sea through CN1 (Fig. 6). Despite no dating age in Pit I is available for this 
sediment layer, it should be earlier than the OSL dating age 1500 (+100) AD at depth 37 cm, which 
coincides with the medieval flooding. Besides, at the depth 62—80 cm in Pit II, the radiocarbon age 
of the plant oddments (871 (+45) a BP at 70 cm, and 800 (+40) a BP at 75 cm), also refer to this 
flooding period considering the bias. Perhaps Pit II was far from the river mouth at that time, as a 
result, the grain size stays low in this layer. It can be speculated that the Lop Sea was not in a stable 
lacustrine deposit situation according to the fact that the grain size of the upper and lower layer in 
Pit I is small, and no plant debris with grain size less than 80 cm and more than 62 cm is found in 
Pit II. 


5.5 Evolution episodes of the drainage system in the eastern Tarim Basin 


From the discussion above, six episodes in the evolution of the drainage system in the eastern Tarim 
Basin during the historical period could be deduced (Fig. 9). 
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Fig. 9 Sketch map of the inferred evolution episodes of the lower Tarim River and Lop Nor over the past 2000 a. 
The water area of the lower Tarim River was labeled in each subfigure. The cylindrical equal-area projection was 
used to project the maps. The water area is in blue polygons, and their outlines are dotted lines if the location was 
speculated without direct field evidence. (a), 0-400 AD was based on topographic data and the locations of the 
ruins; (b), 1260-1450 AD was based on the discussion above. (c), 1760 AD was adopted from the Qianlong 
Thirteen-Row Atlas; (d), 1879-1901 AD was integrated from the expedition reports of Hedin (1940); (e), 1906 AD 
was recovered by the descriptions of Huntington (1907a); and (f), the map was the surveying results of Hörner and 
Chen (1935). 


In the Loulan period, there might be a large brackish lake. Despite the determination of the 
existing time for this 790-m Lop Sea relies upon more additional dating materials, it could be 
suggested that the area (Fig. 9a) was the maximum possible space in the Loulan period. Moreover, 
there should have been rivers flowing by the ruins of LA, LE, and LK (Fig. 6). Under the reduced 
humidity and heat of the westerlies in central Asia during the late Holocene (Chen et al., 2008), or 
anthropogenic activities (Mischke et al., 2017), the volume of the Tarim River was insufficient to 
maintain its previous length. Thus, Lop Sea Basin might be desiccated with the extinction of the 
Loulan civilization. During a flooding period (1260-1450 AD), the river recovered its length 
flowing eastwards, carrying the detrital components in the Loulan depression and discharging into 
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the Lop Sea through CN1 (Fig. 9b), where the water area of the Lop Sea was derived by the contour 
through Pit II. After the medieval flooding, at least before 1760 AD, the Tarim River reduced its 
length again and terminated at Qing Lop Nor (Fig. 9c). Before 1879 AD when Przhevalsky first 
investigated the Tarim Basin, the Tarim river had changed its channel and flowed southward into 
Kara-buran and Kara-koshun (Fig. 9d). With the formation of CN2 between 1901 and 1906 AD, 
the Lop Sea got water supply again (Fig. 9e). Later in 1931, the Tarim River diverted at Temenpu 
and flowed through Kuruk-daria eastward into Lop Sea Basin (Fig. 9f) (Hörner and Chen, 1935; 
Chen, 1936; Hedin, 1940). 

As we know, the diversion of drainage channels was general in arid central Asia. A similar 
diversion was documented by long time series of remote sensing images of Cherchen-daria, one of 
the tributaries of Kara-buran (Tetima) at the southeastern Tarim Basin (Abdimijit et al., 2014). 
However, more researches are needed to identify the details and mechanisms between these 
evolution episodes. 


6 Conclusions 


Based on the current results, we may conclude that the Lop Sea was not on a permanent lacustrine 
deposit condition over the last 2000 a. The terminal lake of the Tarim River, Lop Nor, changed its 
location and area according to the water supply. Interpretations of the sediments in the Lop Sea 
must be treated with care. 

The lake level could hardly exceed 790.0 m in the Loulan period, thus the corresponding 
maximum lake area was about 8000.0 km. The evolution between the Loulan period and medieval 
flooding period was not yet clear, but the Lop Sea might be desiccated with the demise of the 
Loulan civilization. After the medieval flooding, Lop Sea did not receive drainage of the Tarim 
River until 1901-1906 AD. Lop Nor in Qing Dynasty was in a region around Tikenlik. The water 
area of the lower Tarim River was decreasing in the last 2000 a. 

Further fieldwork would be needed to figure out the paleo environment and human activities in 
the Lop region, and the DEM provides a reference for the sample site in the Lop Sea. 


Acknowledgements 


This work was supported by the National Natural Science Foundation of China (42071313, 41671353, U1303285, 
41571363, 41431174, 61471358, 41201346, 41301394, 41301464). 


References 


Abdimijit A, Alimujiang K, Alishir K, et al. 2014. Monitoring the water area changes in Tetima-Kanglayka Lakes region over the 
past four decades by remotely sensed data. Journal of Lake Sciences, 26(1): 46—54. (in Chinese) 

Abdimijit A, Alimujiang K, Alishir K, et al. 2016. Evolution of small lakes in lower reaches of Tarim River based on multi-source 
spatial data. Geographical Research, 35(11): 2071—2090. (in Chinese) 

Chen F, Yu Z, Yang M, et al. 2008. Holocene moisture evolution in arid central Asia and its out-of-phase relationship with Asian 
monsoon history. Quaternary Science Reviews, 27(3—4): 351-364. 

Chen P C. 1936. Lop Nor and Lop Desert. Acta Geographica Sinica, 1: 19—49. (in Chinese) 

Coianiz L, Ariztegui D, Piovano E L, et al. 2015. Environmental change in subtropical South America for the last two millennia 
as shown by lacustrine pigments. Journal of Paleolimnology, 53: 233-250. 

Dong Z, Lv P, Qian G, et al. 2012. Research progress in China's Lop Nur. Earth Science Reviews, 111(1—2): 142-153. 

Fritz T, Balss U, Bamler R, et al. 2012. Phase unwrapping correction with dual-baseline data for the TanDEM-X mission. 
Geoscience and Remote Sensing Symposium, 5566-5569. 

Gonzalez J H, Bachmann M, Scheiber R, et al. 2010. Definition of ICESat selection criteria for their use as height references for 
TanDEM-X. IEEE Transactions on Geoscience and Remote Sensing, 48(6): 2750-2757. 

Guo S, Li D. 2003. A survey of Felix da Rocha's scientific work in China. Palace Museum Journal, 4: 005. (in Chinese) 

Guo S, Li D. 2005. P. José de Espinha, imperial first supervisor in charge of the Western Affairs and Astronomy. Journal of Inner 
Mongolia Normal University (Philosophy and Social Science), 2: 008. (in Chinese) 


ZHANG Tingting et al.: A study on historical location and evolution of Lop Nor in China... 


Han C, Lv G. 2006. The changes of Lopliks living in the east of the Tarim Basin and their environment since the Qing Dynasty. 
Journal of Chinese Historical Geography, 21(2): 60—66. (in Chinese) 

Hedin S. 1905. Scientific results of a journey in central Asia, 1899-1902 Vol. 2 Lop Nur. Lithographic institute of the General 
staff of the Swedish army, 257-369. 

Hedin S. 1940. The Wandering Lake. New York: E P Dutton and Company Inc., 80—112. 

Hörner N G, Chen PC. 1935. Alternating lakes. Some river changes and lake displacements in central Asia. Geografiska Annaler, 
17: 145-166. 

Huntington E. 1907a. Lop-Nor. A Chinese lake. Part I. The unexplored salt desert of Lop. Bulletin of the American Geographical 
Society, 39: 65-77. 

Huntington E. 1907b. Lop-Nor. A Chinese lake. Part II. The historic lake (Lop-Nor). Bulletin of the American Geographical 
Society, 39(3): 137-146. 

Ivory S J, Russell J. 2018. Lowland forest collapse and early human impacts at the end of the African Humid Period at Lake 
Edward, equatorial East Africa. Quaternary Research, 89(1): 7-20. 

Jia H, Liu J, Qin X. 2011. Early Holocene climate changes and agricultural activities inferred from spore-pollen of Lop Nur. 
Journal of Jilin University (Science Edition), 41: 181—194. (in Chinese) 

Jiang H, Shevenell A, Yu S, et al. 2015. Decadal-to centennial-scale East Asian summer monsoon variability during the Medieval 
Climate Anomaly reconstructed from an eastern Tibet lacustrine sequence. Journal of Paleolimnology, 54: 205-222. 

Li B, Ma L, Jiang P, et al. 2008. High precision topographic data on Lop Nor basin's lake "Great Ear" and the timing of its 
becoming a dry salt lake. Chinese Science Bulletin, 53: 905-914. 

LiG, Tang X, Zhang C, et al. 2017. Multi-criteria constraint algorithm for selecting ICESat/GLAS data as elevation control points. 
Journal of Remote Sensing, 21: 96—104. (in Chinese) 

Li K, Qin X, Zhang L, et al. 2018. Hydrological change and human activity during Yuan—Ming Dynasties in the Loulan area, 
northwestern China. Holocene, 28(8): 1266-1275. 

Liu C, Zhang J F, Jiao P, et al. 2016. The Holocene history of Lop Nur and its palaeoclimate implications. Quaternary Science 
Reviews, 148: 163-175. 

Lu Z. 2008. Folk maps of Qing Dynasty. Land and Resources, 12: 58—61. (in Chinese) 

Lucieer A, de Jong S M, Turner D. 2014a. Mapping landslide displacements using Structure from Motion (SfM) and image 
correlation of multi-temporal UAV photography. Progress in Physical Geography, 38(1): 97-116. 

Lucieer A, Turner D, King D H, et al. 2014b. Using an Unmanned Aerial Vehicle (UAV) to capture micro-topography of Antarctic 
moss beds. International Journal of Applied Earth Observation and Geoinformation, 27: 53—62. 

Lv H, Xia X, Liu J. 2010. A preliminary study of chronology for a newly-discovered ancient city and five archaeological sites in 
Lop Nor. Chinese Scientific Bulletin, 55: 237—245. (in Chinese) 

Ma C M, Wang F B, Cao Q Y, et al. 2008. Climate and environment reconstruction during the Medieval Warm Period in Lop Nur 
of Xinjiang, China. Chinese Science Bulletin, 53: 3016-3027. 

Mancini F, Dubbini M, Gattelli M, et al. 2013. Using unmanned aerial vehicles (UAV) for high-resolution reconstruction of 
topography: The structure from motion approach on coastal environments. Remote Sensing, 5(12): 6880-6898. 

Mischke S, Liu C, Zhang J, et al. 2017. The world's earliest Aral-Sea type disaster: the decline of the Loulan Kingdom in the 
Tarim Basin. Scientific Reports, 7: 43102, doi: 10.1038/srep43102. 

Mischke S, Liu C, Zhang J. 2020. Lop Nur in NW China: Its natural state, and a long history of human impact. In: Kalweit S I. 
Large Asian Lakes in Changing World. Berlin: Springer, 207-233. 

Pavlis N K, Holmes S A, Kenyon S C, et al. 2008. An earth gravitational model to degree 2160: EGM2008. EGU General 
Assembly, 10: 13-18. 

Pavlis N K, Holmes S A, Kenyon S C, et al. 2012. The development and evaluation of the Earth Gravitational Model 2008 
(EGM2008). Journal of Geophysical Research: Solid Earth, 117, doi: 10.1029/2011JB008916. 

Pjevtsoff M V, Tong Y, Tong S trans. 2013. The Expedition of Pjevtsoff. Urumqi: Xinjiang People Press. (in Chinese) 

Przhevalsky N, Huang J (trans.). 1999. Moving toward Lop Nor. Urumqi: Xinjiang People Press. (in Chinese) 

Reinhardt C, Wiinnemann B, Krivonogov S K. 2008. Geomorphological evidence for the Late Holocene evolution and the 
Holocene lake level maximum of the Aral Sea. Geomorphology, 93(3—4): 302-315. 

Shao Y, Gong H, Liu C. 2011. Primary interpretation on shorelines of vanished Lop Nur Lake using multi-source SAR data. 
Journal of Remote Sensing, 15: 645—647. (in Chinese) 

Shao Y, Gong H, Gao Z, et al. 2012. SAR data for subsurface saline lacustrine deposits detection and primary interpretation on 
the evolution of the vanished Lop Nur Lake. Canadian Journal of Remote Sensing, 38(3): 267-280. 

Stein A. 1921. Serindia: Detailed Report of Explorations in Central Asia and Westernmost China. Oxford: Clarendon Press, 50- 


JOURNAL OF ARID LAND 


56. 

Sun Y, Jiang L, Liu L, et al. 2016. Generating and evaluating digital terrain model with TanDEM-X bistatic SAR interferometry. 
Geomatics and Information Science of Wuhan University, 41: 100-105. (in Chinese) 

Wallace L, Lucieer A, Malenovsky Z, et al. 2016. Assessment of forest structure using two UAV techniques: a comparison of 
airborne laser scanning and structure from motion (SfM) point clouds. Forests, 7(3): 62. 

Wang Q, Liu R. 2007. Three Surveyed Atlas of Qing Dynasty. Beijing: Foreign Languages Press. (in Chinese) 

Xi G. 1992. Migration history of Lop Nor and new findings. The Western Regions Studies, 4: 5—16. (in Chinese) 

Xia X, Wang F, Zhao Y. 2007. Lop Nur in China. Beijing: Science Press, 1—27. (in Chinese) 

Xu B, Gu Z, Qin X, et al. 2017. Radiocarbon dating the ancient city of Loulan. Radiocarbon, 59(4): 1215-1226. 

Xu H, Ye M, Song Y. 2005. The dynamic variation of water resources and its tendency in the Tarim River Basin. Acta Geographica 
Sinica, 15: 467—474. (in Chinese) 

Yang X, Liu Z, Zhang F, et al. 2006. Hydrological changes and land degradation in the southern and eastern Tarim Basin, Xinjiang, 
China. Land Degradation and Development, 17(4): 381-392. 

Yuan G, Yuan L. 1998. An approach to the envionmental changes in Lop-Nur history. Acta Geographica Sinica, 53: 80-89. (in 
Chinese) 

Zhang J, Wang J, Ma H, et al. 2008. Streamflow variations of four source streams to mainstream of Tarim River, Xinjiang and 
water transportation in 2006. Journal of Glaciology and Geocryology, 30: 569-577. (in Chinese) 

Zhang J F, Liu C L, Wu X H, et al. 2012. Optically stimulated luminescence and radiocarbon dating of sediments from Lop Nur 
(Lop Nor), China. Quaternary Geochronology, 10: 150-155. 

Zhou Y. 1978. Discussion on the issue of Lop Nor migration. Journal of Beijing Normal University (Natural Science), 3: 34—40. 
(in Chinese) 


